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ABSTRACT 


Induced aerodynamics from thrust vectoring are investigated in the present 
study by a computational fluid-dynamic method. A thin-layer Reynolds-averaged 
Navier-Stokes code with multiblock capability is used. Jet properties are specified 
on the nozzle exit plane to simulate the jet momentum. Results for a rectangular 
jet in a cross flow are compared with data to verify the code. Further verification 
of the calculation is made by comparing the numerical results with transonic data 
for a wing-body combination. Additional calculations are performed to elucidate 
the following thrust vectoring effects: 

• thrust-vectoring effect on shock and expansion waves, 

• induced effects on nearby surfaces, and 

• thrust-vectoring effect on the leading-edge vortex. 

These calculations are performed on a body with a rectangular nozzle and a 
thick delta wing of 63-deg sweep. The results show that a vectored jet induces 
flow acceleration and lower pressure in the region with flow expansion. Flow 
deceleration and higher pressure exist in the region of flow compression. When a 
jet is vectored downwards, incremental lift, drag and negative pitching moment 
are induced. Furthermore, induced flow acceleration from a vectored jet can 
change the strength and location of the shock and expansion waves. These 
results depend significantly on the magnitude of jet static pressure. Calculation 
with the delta wing indicates that vectored thrust makes leading-edge vortices less 
diffusive. 
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1. INTRODUCTION 


Thrust vectoring has been considered as a means to enhance 
maneuverability and aerodynamic performance of a tactical aircraft (Refs. 1-17). 
This concept usually involves the installation of multifunctional nozzles. With the 
multifunctional nozzles, the engine thrust can be changed in direction without 
changing the attitude of the aircraft. A change in the thrust direction can induce 
extra aerodynamic forces on the aircraft. Therefore, this device can be used for 
lift-augmentation as well as stability and control purposes. When the jet stream 
is deflected in the longitudinal direction, the lift force and pitching stability can be 
manipulated, while yawing stability can be controlled by directing the thrust in 
the lateral direction. 

The augmented lift force from thrust vectoring can also be used for short 
takeoff and landing (STOL) purposes (Ref. 17). This will allow highways, small 
airports, and bombed runways to be bases for fighter operations (Ref. 10), not to 
mention the safer landing operations on aircraft carriers provided by using this 
concept. Also, the induced aerodynamic forces from thrust vectoring devices allow 
the aircraft to operate in the high angle-of-attack flight regime (Ref. 14). Most of 
the current fighter airplanes are limited in high-angle-of-attack maneuverability 
because of degraded control power. The primary reason for the control power 
deterioration is that massive flow separation occurs at high angles of attack on the 
lee side of the airplane. Utilizing the propulsion system and vectoring the 
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exhaust, powered control can provide auxiliary high-angle-of-attack maneuvering 
capability (Ref. 6). 

On the other hand, the pitching and yawing moments generated from these 
kinds of devices can be used to enhance or replace the elevator and rudder control 
power. Test results show that afterbodies of various aircraft comprised 20% to 
25% of the total model length but produced 38% to 50% of the total aircraft drag 
(Ref. 3). Up to half of the afterbody drag results from adverse interference and 
pressure drag in the afterbody region. One way to eliminate the large adverse 
interference effect is to reduce the size of, or simply eliminate, the horizontal and 
vertical tails. This would require some other methods of achieving the necessary 
control power (Ref. 18). A well-designed thrust- vectoring nozzle providing the 
required control power can replace the roles of horizontal and vertical tails. 
Without the adverse interference among control surfaces and afterbody, the total 
drag of the aircraft can be reduced. 

In an attempt to provide vectoring capability by deflecting the jet stream, 
several devices have been considered. Some of the proposed devices are internal 
turning vanes, external turning flaps, gimbal nozzles, and ramp expansion nozzles. 
These devices produce the desired aerodynamic forces on the aircraft by directing 
the thrust in the longitudinal and/or lateral directions. To determine these forces, 
experimental investigation has been under way for decades (Refs. 2, 3, 5, 12, 14, 
19). To effectively incorporate this technology into future aircraft, the propulsion 
system must be efficiently integrated into the airframe. This involves determining 
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the characteristics of afterbody flows in the presence of a jet exhaust. Advanced 
computational methods are becoming increasingly useful in providing analyses of 
this nature. Tasks associated with thrust vectoring computation include induced 
surface pressure calculation, jet deformation effect, crossflow velocity distribution, 
and induced shock wave at high speed. To fully understand the flow properties of 
thrust vectoring, interaction between the deflected jet and external flow should be 
studied first. 

Experiments show that a crossflow jet injected from a surface produces low 
pressures on the surface downstream of the jet exit. There is also a high pressure 
region due to crossflow deceleration (Ref. 20). The trajectory of the jet mainly 
depends on the jet-to-freestream speed ratio, the nozzle exit geometry, the jet 
injection angle, and the jet exit profile (Refs. 21-23). Experiments also show that 
large-scale rolled-up vortices are generated in a crossflow jet (Ref. 24). These 
phenomena are expected to be present in the thrust vectoring operation. 

Reported numerical simulations of crossflow applications are limited to 
simple configurations (Refs. 25-31). For configurations having multifunctional 
nozzles, the flow field about the nozzle exit becomes very complicated due to the 
interaction among the turning devices, the freestream, and the vectored jet. Jet 
mixing as well as shock waves induced by the jet turning will be present. For 
high nozzle pressure ratio (NPR) or supersonic nozzles, the shock waves inside the 
jet stream will affect the jet deformation. Also, pressure differences across the jet 
can produce vortices, which will induce jet roll-up. In addition, vectoring the 
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thrust may create flow expansion in the external region around the nozzle exit. 
Flow expansion in this region induce shock waves at a transonic speed. With all 
the complicated flow physics induced from the vectored thrust, methods beyond 
crossflow prediction are needed for thrust vectoring computation. As a general 
goal of this study, numerical methods resolving induced aerodynamic effects-- 
namely, pressure distribution change on the control surfaces, shock location 
variation around the jet exhaust region, and the effect on vortices-should be 
developed. 

One thrust vectoring computational method based on a panel method was 
developed in Reference 32. Due to the limitation of flow singularity methods, 
important flow characteristics such as shock waves and shock cells inside the jet 
stream can not be predicted. Another drawback of panel methods is the absence 
of the energy (or pressure) term in the computation. Therefore, the jet stream is 
frequently modeled as a plume or wake, requiring empirical methods to determine 
the jet shape. To fully understand the induced aerodynamic effects from the jet 
stream deflection, methods other than a panel approach would be needed. 

Methods in Computational Fluid Dynamics (CFD) are more suitable for the 
present problem. In more recent simulations of unvectored jets, flow 
characteristics at the nozzle exit plane were specified and the exhaust plume was 
calculated based on a Navier-Stokes method (Ref. 33). However, calculation for an 
airplane configuration with a vectored jet has not been reported in the literature. 
Usually flow properties involved in CFD are density, velocity, energy, and/or 
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pressure. Boundary conditions can be specified in terms of these quantities. For 
example, on solid boundaries, velocity components can be specified as zero to 
simulate the no-slip boundary condition. For thrust vectoring computation, jet 
exit conditions can be specified on the nozzle exit plane. The jet deformation or 
trajectory can be calculated instead of using empirical formulas to represent the 
jet effect. 

In the present study of thrust vectoring aerodynamics, nonaxisymmetric 
nozzle with turning vanes are investigated with a multiple-block thin-layer 
Reynolds-averaged Navier-Stokes code, CFL3D (Ref. 35). The numerical grid 
generation program EAGLE (Ref. 34), with the options of algebraic and elliptical 
grid generation systems, has been used to produce grids in this study. To properly 
simulate the jet interaction effect, the calculation is carried out inside the jet 
region as well as the flow field around the nozzle and the airplane configuration. 
This approach captures the jet deflection trajectory through the thin-layer 
Reynolds averaged Navier-Stokes equations. Shock cells inside the jet stream and 
the interaction between the jet and the freestream can also be predicted. For the 
nozzle of a turning-vane configuration, the vanes are not directly modeled in the 
present study. Instead, the turning angle, velocity components, and pressure of 
the jet are specified as boundary conditions on the exit plane. Under most 
circ um stances, available test conditions inside the thrust vectoring nozzles are 
total temperature and total pressure. To convert total pressure into CFD 
variables, isentropic expansion formulas are used to determine the jet exit 
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pressure. The temperature effect is not considered at present. The Navier-Stokes 
approach of CFL3D (Ref. 35) is based on a finite volume scheme and will be used 
here for a thrust vectoring nozzle configuration. The finite volume scheme 
features momentum conservation in the computation. Therefore, this scheme is 
suitable to calculate the jet momentum transfer. 

To verify the jet flow calculation, a 90-degree crossflow jet is treated first. 
The pressure distribution is compared with experimental measurements. A 
generic airplane with a rectangular nozzle is used to demonstrate the calculation 
of the induced pressure and shock location changes due to a vectored jet. 
Computational results of this configuration with various NPR’s and deflection 
angles will be presented and compared with test results. 

To investigate the vectored-jet effect on the leading edge vortex, a 63-degree 
delta wing integrated with trapezoidal nozzles is used. Pressure distributions and 
leading-edge vortex filaments under the influences of a straight jet and a vectored 
jet will be presented. Vectored-jet influence on vortex breakdown will also be 
discussed. 
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2. THEORETICAL APPROACH 

Generally, tasks in developing CFD programs include derivation of the 
governing equations, solution of the flow equations, accuracy and stability analysis 
of the scheme, grid generation, etc. In this section, outline of the theories as used 
in the CFL3D code will be presented. The grid generation method used in this 
study will also be discussed. 

The program CFL3D utilizes thin-layer Reynolds averaged Navier-Stokes 
equations, enabling detailed flow properties to be obtained. In solving the 
governing equations, the program is based on a finite volume approach. In 
addition, the accuracy of the solution depends on the smoothness of the grids. 

In spatially discretizing equations into finite difference forms, central 
difference is usually used. This formulation ensures a second order accuracy. 

From the consideration of accuracy, a central difference is better than either the 
backward or forward difference. However, in some particular situations, backward 
difference has to be used to satisfy the physics of the flow field. For example, in a 
supersonic flow region, the Mach lines drawn from a point divide the flow into 
zones of influence. Properties at this point can only be affected by flow in the zone 
of influence. Therefore, while formulating differentiation in difference forms, only 
the points in the zone of influence should be taken into consideration. In this 
situation, the backward difference scheme has to be used. For subsonic flow, 
characteristics at each point can be affected by any other points in the flow field. 
Thus, the central difference scheme is suitable for subsonic flow calculation. 
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In transonic flow computations, the formulation of the difference equations 
becomes more difficult due to combination of subsonic and supersonic regions. As 
stated earlier, in the subsonic region, the central difference scheme should be 
used. On the other hand, in the supersonic region, backward difference has to be 
applied. Upwind difference was developed to solve this problem. Take the 
differentiation dF/dE, as an example. The formulation can be written as 


dF _ ir F i+l ~ F i-1 

d$ ^ 2A^ 


•) + (l- pX 


Fi - Fj.i, 


( 1 ) 


p = 1 subsonic 

p = 0 supersonic 


Therefore, central difference is used in the subsonic region. In a supersonic 
region, backward difference is used. 

In CFL3D, options of central difference with added dissipative terms, flux 
difference splitting, and flux vector splitting are provided to separate the flow field 
into regions with different flow characteristics. Details of upwind differencing 
used in the program will be discussed later. In addition, CFL3D has the 
capability of multiblock computation. The development of multiblock computation 
enables CFD to be applied to complicated configurations. Also, this concept makes 
the grid generation an easier task. For example, an individual block can be 
formed in the region where smooth continuation of the grids is not possible. 
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Therefore, smooth grids can be constructed in each block. Flow properties can 
also be solved accurately in each block. Transmission of flow characteristics 
among blocks can be achieved through interpolations. 

2.1 Governing Equations 

From the conservation of mass, momentum, and energy, the Navier-Stokes 
equations can be written as 

^ + JL (F - F v ) + — (G - G v ) + A (H - H v ) = 0 (2) 

dt dx dy dz 


where 


Q = 


\ 

P 

pu 

pv 

pw 


Flux in the x-direction: 


/ A 


(o ) 

pu 

2 


T-xx 

p ^ + P 

puv 

<i 

ii 

'txy 

puw 


T xz 

v (e + p)u y 


^ u ^xx + v ^xy + wi; xz ” 


(3) 
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Flux in the y-direction: 


G = 


pv 

puv 

pu 2 + p 
puw 
(e + p)v 


G„ = 


l xy 

c yy 


xy 

ux 




yx + VT yy + wX yz 


^y 


( 4 ) 


Flux in the z-direction: 


H = 





'o ^ 

pw 



pwu 



X zx 

pwv 


H v = 

x zy 

pw 2 

+ P 


x zz 

^(e + 

P)w y 


^ X zx + vx zy + ^ x zz ~~ 4z^ 


( 5 ) 


where e is the total energy per unit volume defined as 

e = 5 — + -ip(u 2 + v 2 + w 2 ) 

y - 1 2 


The shear stress x x . x ^ and heat flux q x . are defined as 


M Q 


X X;X' = 


,A J R, 


p( 


t s„ 

3x; 9xj dxfc 


a -Ik?! 
qxi P r 9x 
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2.2 Coordinate Transformation 


Because of the irregular grid spacing in the physical domain (x-, 
y-, z-coordinates), the order of computational accuracy is difficult to predict. From 
a scheme-stability point of view, it is better to perform the calculation in a 
generalized coordinate domain (so-called computational domain). The generalized 
coordinate transformation can be described as 

% = £(x, y, z, t) 
ti = q(x, y, x, t) 

C = C( x > y> z > t) 


The Jacobian of this generalized transformation is obtained as 


J = D. 0 

d(x, y, z) 

= [x^y^ - y^) - y^(x ri z T - x^) + z^(x n y^ - y^)]" 1 


The governing equation becomes 


+ JL (F - F v ) + JL (G - G v ) + JL (H - H v ) = 0 
at an ac 


( 8 ) 


where 
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Q = Q/J 


{• - F v = Mi [4 X (F - F v ) <■ ? y (G - G v ) ♦ 4 Z (H - H v ) * ?Q] 
j 

G - G v = iXlL [i\ X (F - F v ) + f\ y (G - G v ) + t[ Z (H - H v ) + f\ t Q] 
J 

H - H v = JZii K„(F - F v ) * Cy(G - G v ) t £ z (h - H v ) * £ t Q] 
j 

(K x , K y , K z , Kt) = (K x> K y , K z , K t V|VK| 

|VK| =. (K^ ♦ Ky <■ K\) m 

S = (^ x u + ^ y v + E, Z w)/|VK| 
v = (tj x u + T) y v + T| z w)/|VK| 
w = (Cx u + CyV + C Z W V| VK[ 

K = (£, T|, CX respectively 

The derivation of Equation (9) can be found in Appendix A. 
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2.3 Thin-Layer Navier-Stokes Approximation 

It is known that the velocity and pressure vary rapidly near a solid surface 
in the direction normal to the surface. On the other hand, flow characteristics are 
generally smooth along the surface, unless singularities exist. Retaining only the 
viscous gradient terms corresponding to the direction normal to the surface can 
save computing resources, yet maintain enough computational accuracy. In this 
manner, the thin-layer Navier-Stokes equations can be obtained as 

<!§ ♦ A (?) + JL (G) + A (ft - H v ) = 0 <10) 

3t 31; Sn 3C 

Certainly, viscous gradient terms can also be retained in the and T| directions 
because, in CFL3D, solid surfaces can be placed on any block surfaces. The 
laminar viscous terms can be included in all tj, £ directions simultaneously. 
However, turbulence viscous terms can only be added at most in two directions 
simultaneously. 

It should be noted that the turbulent equations are derived with a time- 
averaging procedure. By decomposing every flow variable into the sum of average 

value and fluctuation, such as Uj = Ui + u[ , time-averaging of momentum 
equations would produce unknown terms such as pu[uj . These new unknowns 
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(called Reynolds stresses) need to be specified to close the system of equations in a 
technique called turbulence modeling. In the present application, turbulence 
modeling is accomplished by an algebraic method in which these Reynolds stresses 
are specified in a form similar to the viscous terms. As a result, the final form of 
the equations to be solved is the same as Eq. (10) except that the coefficient of 
viscosity (p) is to be replaced by p + e, where e is the eddy viscosity. In Section 
2.8, the way e is calculated will be summarized. 

2.4 Approximate Factorization 

The thin-layer Reynolds-averaged Navier-Stokes equations are solved by 
three-factor approximate factorizations. For example, at time step n + 1, Equation 
(10) can be written as 


[_L + 8p — 

JAt q 3Q 



3G 

3Q 



^I)] n AQ = -R n 


(ID 


where R n is the residual. Equation (11) can be factorized into three equations 
along r|, £ sweep directions as follows: 


[J_ + 8p —] n [l + JAtS 
JAt ^ 3Q 


f!!±] n [I + JAt5r(M 

3Q ^ 3Q 


^5l)] n = -R n 
3Q 


[_L + ^] n AQ* = -R n 

JAt ^ 3Q 


( 12 ) 
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( 13 ) 


[_L + 6. <*ij n AQ 
JAt 3Q 


** = (__ )AQ* 
JAt 


JAt 


5 C ( 


3H 


^I)] n AQ = (-L)AQ 
9Q 9Q JAt 


* ★ 


(14) 


The solution is then updated as 


Q n+1 = Q n + AQ 


(15) 


2.5 Upwind Differencing Scheme 

Three distinct upwind difference schemes are provided in this program; e.g., 
flux difference splitting, flux vector splitting, and central difference with explicit 
artificial damping. The central difference scheme used in this program is the 
same as that used in Reference 36 by adding the negative of the following explicit 
dissipative term to the governing equations: 

Dqijk = D x qjjk +D y qjjk + D z qjjfc 

D x<lijk = d i+l/2j,k “ d i-l/2j,k 

Dyqjjk = d ij+l/2,k " d ij-V2,k (17) 

DzOijk = d ij,k+l/2 " d ij,k-l/2 
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( 18 ) 


d i+l/2,i,k = 1+ ^~ [ £ i?i/2j,k^- t -lj> k qi J. k ^ 

At 


“ e i + 1/2 j ,k^ q i + 2 j ,k 3 fli+lj,k + 3( lijk fli-lj,k] 


Define 


v » k _ Pj+lj,k 2 ^ijk + Pj-lj.k 
|Pi+lj,k + 2Pijk + Pi-lj,k| 

e i+l/2,j,k = k (2) max(v i+ i )jik , v ijk ) 


(19) 


( 20 ) 


e 


(4) 

i+172,j,k 


max[0,(k (4) - zZ } V2,i,k )] 


( 21 ) 


where k (2) and k (4) are artificial damping coefficients. Typically, values of k (2) 
and k*- 4 -* are 1/4 and 1/256, respectively. 

In Equation (20), the artificial damping sensor Ujj k is used to sense the 
existence of shock waves. Theoretically, the pressure distribution is smooth in the 
flow field if no shock waves exist. Therefore, u- k is close to zero. The central 
difference is used in this situation. If a shock occurs in the flow field, the pressure 
rises across the shock. The sensor has a value close to one; thus, backward 
difference becomes the dominant term. 
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The schemes of flux difference and flux vector splittings are summarized 

below. 

2.5.1 Flux Difference Splitting: The generalized fluxes F, G, and H are split 
into forward and backward contributions according to the signs of eigenvalues of 
the Jacobian matrices. Eigenvalues of the Jacobian matrices determine the 
direction of the propagation waves. If an eigenvalue of the matrix is positive, the 
wave corresponding to this eigenvalue is moving toward the positive direction. 
Therefore, backward difference should be used. Otherwise, if the eigenvalue is 
negative, forward difference should be used. For example, consider the flux in the 
£ direction, 

A = — = TAT -1 = T(A + + A")T _1 (22) 

3Q 

where A* = - ± 1M . A are diagonal matrices formed from the eigenvalues of A. 
2 


The left-hand side of Equation (13) can be expressed as 

[ 1 + 5 ? i?.] n AQ* - T[_i_ + 5:a + + 8tA']T _1 AQ* (23) 

JAt ^ 3Q JAt ^ * 

In the same manner, the flux in the t\- and ^-directions can be split accordingly. 
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2.5.2 Flux Vector Splitting: The generalized fluxes F, G, and H are split 

into forward and backward contributions according to the ratio of flow speed and 
speed of sound. 

For example, consider the flux in the ^-direction in Equation (10): 


8^F = 5jTF + + 5^F' 


(24) 


where 5^ and 5^ denote the backward and forward differences. 


Define: 


= u/a and u = u/|V^| 


(25) 


for supersonic flow (|M^| > 1) 


F + = F, F“ = 0, > 1 

F + = 0, F“ = F, < -1 


and for subsonic flow ( |M^ | < 1) 
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F- = 


_ m\ 


r 

mass 

± 2a)/Y + u] 

Cass^" 5 ± 2a)/ ? + v] 
Cass^zC-u ± 2a)/y + w] 
f* 

energy 


C» - d 2/4 < m e, - f > 

^energy " * 2 <?-» 53 + 2a 2 ]/( Y 2 + l) 

+ (u 2 +v 2 +w 2 )/2} 

k x = ^ x /|V^|, k y = ^y/|V^|, k z = ^ Z /|V^ | 


( 27 ) 


The flux in the r\- and ^-directions can be obtained in the same manner. 
Derivation of the flux vector splitting formulation can be found in Appendix B. 
The split-flux difference is implemented as a flux balance across a cell as 


8^F = 5^F + + 5£ F~ 


( 28 ) 


Since = Arj = A£ = 1 , the split-flux difference can be written as 
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Q is evaluated at cell interface by upwind-biased interpolation. 


Q 1 = + — <t>^[(l - kjOV^ + (1 + k^)A^]Qi 

i+_ 4 

2 


Q + i = Qi + i - — + k^)V^ + (1 - k^)A^]Q i+ i 

i+_ 4 

2 


where 


A^Qi = Qi + i - Qi V^Qi = Qi - Qi_i 

(|)t = 0 for the first-order upwind, and 

<(>^ = 1, k^ = -1 for the second-order upwind, and 


<t>£ = h k E, 


i for the third -order upwind. 
3 


2.6 Multiblock Operation 


( 30 ) 


(31) 


Usually, an attempt to generate a single grid for a complicated geometry 
results in highly skewed grids, which in turn result in inaccurate calculations 


(Ref. 38). For example, the grid system of the unwrapping type can only describe 
a simple configuration, such as wing alone or body alone cases. Applying the 
same grid system to a wing-body configuration, grid smoothing is needed. Also, 
this grid system usually produces coarse grids on the fuselage while using the 
same number of grid points on the wing in the unwrapping direction. For a 
configuration with integrated nozzles, the single-block calculation (gnd generation 
and flow solver) cannot produce accurate results. Therefore, multiblock 
algori thm s become inevitable in the calculation of complex configurations. 

In the multiblock calculation, the computational domain is divided into 
several hexahedral blocks. Basically, the program solves the thin-layer Navier- 
Stokes equations in individual blocks. Boundary conditions on boundary surfaces 
of each block have to be specified. If a boundary surface is designated as a patch 
surface, provision must be made for communication across this block interface. 

Generally, the zonal approach bears a number of advantages. Firstly, it 
avoids the difficulty of generating three-dimensional grids for different types of 
configurations with appropriate clustering to capture high-gradient flow 
parameters (Ref. 38). Secondly, zonal methods allow different types of grid 
topologies to be used to improve mesh efficiency. Thirdly, flow solvers can be 
written in a structured manner; that is, solution algorithm remains unchanged 
with appropriate boundary conditions on the block faces to calculate different 
configurations. 
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For global computation, a mechanism to transfer information between 
neighboring blocks must be constructed. Two kinds of techniques, patched-grids 
(Ref. 39) and overlapping-grids technique (Ref. 40), have been widely used in 
multiblock applications. In a patched-grids approach, blocks meet at a common 
interface. Three types of grid arrangements can be used for patched-grids 
approach (Ref. 39): 

(a) one-to-one matching, in which grid lines match across the interface; 

(b) MxN-to-one matching, in which grid lines of a cell face match the 
outer contour of MxN cell surface; and 

(c) arbitrary matching, in which there is no restriction on the grid lines’ 
matching. 

In the overlapping-grid approach, an overlapping region is required to provide 
matching solutions across boundary interface. The requirement of a common 
interface between blocks places restriction on the grid generation process for the 
patched-grid approach. For the overlapping-grid approach, a great amount of 
interpolation in the overlapping region is usually required. 

CFL3D uses a grid-block patching algorithm based on generalized 
coordinates. For example, consider two blocks with a common interface. Suppose 
the spatial flux into block 1 at the interface is desired. Two of the four data 
points required to construct the split fluxes at the interface come directly from the 
block 1 values in cells adjacent to the interface. The other two data points can be 
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obtained by interpolation of the block 2 data to block 1 ghost cell centers projected 
from block 1. 

On each interface plane, flow properties at each cell center are obtained by 
interpolating flow properties at cell centers of neighboring cell from adjacent block 
patch surface. To relate the flow properties of block 2 to the image cell center of 
block 1, interpolations are needed. The interpolations are accomplished by 
correlating the cell-center properties in block 2 into a polynomial with constant 
coefficients. In general, a nonlinear polynomial will be formed; therefore, the 
Newton method is used to determine the corresponding flow properties of the 
image plate points projected from block 1. 

2.7 Jet Application 

In the program CFL3D, eight standard boundary conditions are available. 

A list of the standard boundary conditions can be found in Appendix C. Special 
purpose boundary conditions can be added following the multiblock coordinate 
orientation. In the current work, a jet block has to be formed to simulate nozzle 
exit conditions. The nozzle exit plane is treated as one boundary surface of the jet 
block. A special routine has been developed to specify the nozzle exit conditions 
and include the jet in the computation. On the other boundary surfaces, boundary 
conditions can be specified as patching or in-/out-flow surfaces. 

The nozzle exit conditions are 
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( 32 ) 



Usually, the nozzle pressure ratio (NPR) is given for a certain test 
condition. The nozzle exit pressure can then be obtained through the assumption 
of isentropic expansion: 

= a ♦ ? - 1 m 2 )^ 71 <33) 

Pj 2 

The thrust vectoring angles are determined by the ratios between velocity 
components u, v, and w. The jet deflection angle in the longitudinal direction, 8, 
equals tan' 1 (w/u); the deflection angle in the lateral direction equals tan' 1 (v/u). 

The jet boundary conditions in the computation are set up as cell center 
boundary conditions. This implies that the boundary conditions are specified at 
the cell centers of the cells adjacent to the boundaries. For example, if the jet is 
moving to the right-hand side, the centers of the cells on the left-hand side of the 
nozzle exit plane are specified with the designated nozzle exit conditions. Usually, 
the jet conditions are assumed to be uniform inside the nozzle. For thrust 
vectoring computation, the flow region starting from the nozzle exit plane can be 
computed according to the Navier-Stokes equations. 
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2.8 Turbulence Model 

The effect of turbulence modeling is accounted for through the concept of 
eddy viscosity. The eddy viscosity is computed by the standard Baldwin-Lomax 
algebraic turbulence model (Ref. 41): 


e 


e«N; w - e sw N; 

X T 2 XT 2 

N sw + N w 


(34) 


where e w and N w are the eddy viscosity and directed distance from the wall. 
In the Baldwin-Lomax two-layer model, the eddy viscosity is given by 


e 


< e inner y ^ ycrossover 
e outer ycrossover < y 


(35) 


where y is the normal distance from the wall and y crossove r sma ^ est va ^ ue 

of y at which values e inner equals £ outer . 

The eddy viscosity of the inner layer is determined by 

e inner = p C 2 1 CO | (36 


where 

C = ky[l - exp (y + /A + )] 


0) | is the magnitude of the vorticity. 
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and 



The eddy viscosity of the outer layer is determined by 


e outer = KC cp F wake F KLEB (y) 


where K is the Clauser constant, and C cp is another coefficient. 
F wake is determined by 


F wake = min 


Ymax F max 

or 


C wk^maxU max 


F max = max(y|co|[l - exp(-y + /A + )] 
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( 37 ) 


The value of 


( 38 ) 
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The value of I^kLEbO^ is 


F KLEB(y) = 
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1 + 5.5 
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CklebY 

ymax 




1-1 


) -> 


(39) 


The constants used in CFL3D are 
A + = 26 
C cp = 1.6 

^KLEB = 03 

c wk = 0.25 

K = 0.4 
k = 0.018 
2.9 Grid Generation 
2.9.1 Surface Grid Generation 

If the boundaries of a surface are specified, the coordinates on the surface 
can be obtained by various interpolation methods. Popular interpolation methods 
include Lagrange interpolation, Hermite interpolation, spline interpolation, etc. In 
the EAGLE program (Ref. 34), the Lagrange interpolation is recommended, 
although other interpolation functions are available. 
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The general form of one-dimensional (£,) Lagrange interpolation is 


N 




N 


Y(£) = E $ n(y)Y n> < t> n^y) = n 

n=l 1 1 «=1 


5 -*• 

U - *1 


( 40 ) 


where N is the number of points in the t, interpolation domain and ^ £ [0, I]. 

For example, for a two-point Lagrange interpolation. 


Y(4) = (1 - i)Yi * |Y2 


(41) 


The general form of two-dimensional (£, ri) Lagrange interpolation is 


N „ M 

y(£, 11) = E ^nCy^^n.ll) + E m) 

n=l 1 m=l J (42) 

N M - 

- EE ^ n^y)Vm(y)Y(^ n» 11 m) 
n=l m=l 1 d 


where 


Vm(-j) = 


M ri - r\ t 

n_ — — 

{=1 11 m “ lit 


M is the number of points in the r| interpolation domain, and r\ e [0, 1]. 
From Equation (42), the coordinates of the grid points can be obtained. 
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2.9.2 Three-Dimensional Grid Generation 

2.9.2.1 Algebraic Grid Generation: In a three-dimensional domain (£,, r\, Q, 
the interpolation becomes 

N j. M 

T(^» T|, 0 = 52 1 1* 0 + 52 'H m» 0 

n=l 1 m=l J 

T r N M , 

+ 52 e t (l)Y(^ n. Ct> - E E d, 0 

t-1 K n.lm-1 I j (43) 

N T , r MT tir 

- E E ♦ n4) 0 n> Tl, Ct) - E E Vm(4)0 t (f>y(^ *1 n» Ct) 

n+1 t=l 1 k m=l t=l J k 

N M T „ r 

* E E E*n<4)¥m4)et(f)ir«n.’lm.Ct) 

n=l m=l t=l 1 d K 

The grids can be produced by using Equation (43). These are the so-called 
algebraic grids. Usually this type of grid generation system produces smooth grid 
coordinates. The main advantages of this algebraic grid generation are simplicity 
and saving of computing time. From Equation (43) the grids are produced based 
on the boundary spacing. Therefore, boundary point spacing is critical in this 
system. For example, in flow calculation the spacing near abruptly changing 
contour should be small. Therefore, the spacing should be small while specifying 
the boundary points for surface generation in the EAGLE. In other words, users 
need good judgement while using the algebraic grid generation system. One way 
to avoid the disadvantages inherent in using the algebraic grid generation system 
is to use the elliptic grid generation system. 
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2.9.2.2 Elliptic Grid Generation: The purposes of elliptic grid generation 
system are controlling coordinate line distribution, and orthogonality in the field. 
Elliptic grids can be generated by solving the Poisson equations: 

V 2 ^ i = P i (i = 1, 2, 3) (44 


in which p 1 is the control function. 

Many types of control functions can be used in Equation (44). Control 
functions to control spacing, orthogonality, curvature, etc., are all different. For 
example, functions to control the spacing can be written as 



i 

Ui 


(45) 


where y is the arc length of a grid line along the boundary surface. 

Detailed formulation of the control functions can be found in Reference 34. 
In solving the Poisson equations, the algebraic grids obtained from the transfinite 
interpolation are used for initial guesses. Then iterative procedures are used to 
solve Equation (44) until satisfactory results are obtained. 
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3. RESULTS 


For the nozzle flow calculation, validation of the computational method will 
be demonstrated first. Subsequently, several aerodynamic phenomena related to 
thrust vectoring will be discussed. They are 

(1) induced pressure variation on the solid surface, 

(2) thrust vectoring effect on expansion and shock waves in the jet, 

(3) jet static pressure effect on the thrust vectoring effect, and 

(4) thrust vectoring effect on the leading edge vortex. 

3.1 Jet Model Validation-a 90-Degree Rectangular Jet 

A physical model simulating the vectored jet flow is first employed in this 
study. Validation of this jet model is necessary before performing the thrust 
vectoring calculation. 

Usually, the process of computational fluid dynamics validation is described 
from two perspectives: numerical- error validation, and physical-model validation 
(Ref. 42). Errors associated with the first aspect are time and space discretization, 
grid refinement and numerical dissipation. To eliminate uncertainty in the 
numerical validation, a simple configuration in steady flow is preferred. The 
second aspect of CFD validation is associated with the governing equations and 
turbulence models. In CFL3D, thin-layer Reynolds-averaged Navier-Stokes 
equations and Baldwin-Lomax turbulence models are used. The Baldwin-Lomax 
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turbulence model is mainly developed for a solid surface such as a wing. This 
model is not appropriate for jet turbulence effect. Therefore, in the first example, 
a laminar-flow case is used. 

An experimental study of a 90-degree rectangular jet in a crossflow was 
presented in Reference 22. The rectangular jet has an aspect ratio of 0.25. The 
jet is injected into the crossflow from a flat plate. A test condition with a jet-to- 
freestream speed ratio (R) = 2.2 is chosen for computation. The jet exit speed in 
the experiment was 225 ft/sec. With a jet-to-freestream ratio of 2.2, the 
freestream speed is about 100 ft/sec. The Reynolds number based on the flat plate 
length is 93600. Because of the low Reynolds number, turbulence effect is not 
considered in the calculation. As mentioned in the conditions of CFD validation, a 
simple-configuration with laminar flow is desired. For this rectangular jet, grid 
skewness is eliminated. 

In this calculation, flow conditions are taken to be 


M M = 

0.1 

a = 0 


U iet = 

0 

V jet = 

0 

W jet = 

= 2.2 M, 

II 

pT 

Pfnf 


Velocity vectors shown in Figure 1 are the calculated velocity components 
slightly away from the flat plate (at cell centers of the cells adjacent to the flat 
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plate). It is seen that the incoming flow slows down as it approaches the jet. The 
jet boundary acts like a solid body placed in the crossflow. This blockage effect 
causes the incoming flow to travel around the jet. The divided external flows 
eventually merge behind the jet region. As a result of flow merging, circulations 
occur in this region. Pressure distribution on the flat plate is shown in Figure 2. 
When the flow approaches the jet, high pressure region is created as a result of 
flow deceleration. When the flow separates and travels around the jet, flow 
accelerates on the sides and downstream of the jet. As a result, low pressure 
regions are created. Pressure distributions from experiment and calculation are 
presented in Figure 3. As shown in this figure, computed results agree reasonably 
well with the measurement. Evidently, this jet model used in the computation is 
capable of simulating the vectored-jet effect. 

3.2 Thrust Vectoring Interference on Wing-Body Configuration 

An investigation was conducted in the Langley 16-foot transonic tunnel to 
determine the induced lift characteristics of the vectored-thrust concept in which a 
rectangular jet-exhaust nozzle was located in the fuselage at the wing trailing 
edge. A sketch showing the external geometry of the model is presented in Figure 
4. The wing has a leading-edge sweep of 50 degrees, streamwise NACA 64A406 
airfoil sections, an aspect ratio of 3.0, and a taper ratio of 0.3. The fuselage has 
rectangular cross sections with rounded comers and had an effective fineness ratio 
of 7.28. The afterbody boattail angle was 12.5 degrees (Ref. 19). As shown in 
Figure 5, thrust vectoring was obtained by using circular-arc turning vanes. Also 
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shown in the same figure, this configuration has a convergent nozzle, which 
implies that the nozzle exit Mach number can not be greater than unity. For a 
certain nozzle pressure ratio (NPR), the nozzle exit pressure can be estimated 
according to the isentropic assumption, providing the jet flow is choked. 

Difficulties associated with modeling the turning vanes are gnd generation, and 
viscous effect related to the rapid geometric changes. In this study, tuming-vane 
mechanism is not modeled in the calculation. In the simulation of the vectored 
thrust, instead of having the actual turning vanes, nozzle exit conditions as 
described in Section 2.7 are used. 

Because of the complexity of the geometry, ten blocks are used to model the 
flow fields around this configuration. Some grid topologies of this configuration 
are shown in Figure 6. As shown in Figure 5, the top and bottom walls occupy 
about 15% of the nozzle height. To account for the wall effects, these walls are 
modeled as solid surfaces. However, the tinning vanes and supporting struts, 
which occupies about 30% of the nozzle exit area, are not included in the 
calculation. Therefore, total exit area in the computation is about 30% larger than 
the actual exit area. With the same jet speed and pressure, the mass flow rate in 
the computation is 30% larger than the actual jet flow rate. A larger jet flow rate 
might cause the vectored jet effects being overpredicted in the computation. 

As shown in Figure 5c, the vane-tuming angle (8^) is used as one of the test 
parameters. From measurement (Ref. 19), effective jet-turning angles depend on 
vane-tuming angles as well as thrust coefficients. Figure 7 shows the effective jet- 
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turning angle under various vane-tuming angles and thrust coefficients. In the 
computation, effective jet-turning angles are used to simulate jet-vectoring effects. 
Jet total pressure is obtained from Figure 8 at corresponding jet conditions. 

Jet conditions used to investigate the jet-induced variation of pressure are 
as follows: 
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Calculated pressure contours on this configuration are shown in Figures 9 and 10. 
The upper surface pressure contours are shown in Figure 9. Comparing these two 
sets of pressure contours, it is seen that the vectored-jet generates lower pressure 
on the upper surface of the afterbody than the straight jet. Furthermore, a lower 
pressure region is observed on the wing when the vectored-jet is incorporated. On 
the other hand, higher pressure regions are generated on the lower surfaces when 
the vectored-jet is on (Figure 10). Computed pressure distributions on the 
afterbody are compared with data in Figure 11. 

On the upper surface, the predicted pressure distribution shows more rapid 
changes around the comer of the ramp (x/X = 0.2 to x/X = 0.5). This is most likely 
affected by not having enough grid points in the boundary layer. When the 
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number of grid points is increased, the flow simulation is improved (Ref. 43). The 
calculation also shows the thrust vectoring effects on the ramp. At the region 
near the exit (x/X = 0.8 to x/X = 1.), a lower pressure distribution is obtained when 
the jet is vectored. This figure also shows that the vectored jet brings the recovery 
shock toward the jet. 

On the lower surface, the recovery shock is pushed forward by the vectored 
jet. A higher pressure distribution is also induced. Sectional pressure 
distributions on the wing are shown in Figure 12. On the wing, higher pressure 
on the lower surface and lower pressure on the upper surface are induced by the 
vectored jet. The recovery shock on the upper surface is also moved backwards by 
the vectored jet. 

As mentioned earlier, computational discrepency occurs at the ramp comer 
region. This region features rapid geometric variation. Expansion waves and flow 
acceleration are expected. One possible reason for the discrepancy is that the 
recovery shock waves are not well captured because of insufficient grid resolution 
in the boundary layer. To correctly predict the shock strength and location, very 
dense grids and/or adaptive grids are needed. The grid density in the current 
calculation is not high enough to resolve the boundary layer flow and capture the 
strong shock waves. A lower Mach number is then used to demonstrate the thrust 
vectoring calculation. This configuration is calculated under the following 
conditions 
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Figure 13 shows the pressure distribution on the center plane of the afterbody. 

The calculated results agree with the experiments well. Again, the calculation 
demonstrates the thrust vectoring effects; i.e., lower pressure on the upper surface 
and higher pressure on the lower surface of the afterbody. Pressure distributions 
on the wing are shown in Figure 14. Again, a more negative pressure distribution 
on the upper surface and a more positive pressure distribution on the lower 
surface are observed. The calculation shows that the vectored-jet induces similar 
thrust vectoring effects on the wing as well. 

3.3 Thrust Vectoring Effect on Shock and Expansion Waves 

After the computational method has been verified, some thrust vectoring 
effects will be elucidated based on computed results. For this purpose, the 
fuselage-alone configuration of Figure 4, as shown in Figure 15, is used to 
calculate the thrust vectoring effect on shock and expansion waves. This 
configuration has rectangular cross sections with rounded comers and has an 
effective fineness ratio of 7.28. This configuration was calculated under the 
following jet conditions. 
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Grid topologies of this configuration are shown in Figure 16. For a boattail 
configuration, flow expansion will occur at the deflecting point. At a transonic 
speed, if the ramp angle is large enough, a supersonic region will be created by the 
flow expansion. As the flow travels downstream, it eventually becomes subsonic. 

A "recovery shock" will occur between the expansion waves and the free stream 
flow downstream. In this section, thrust-vectoring effects on the expansion waves 
and the recovery shock will be discussed. The expansion ramp angle of this 
configuration is 12.5 degrees. With a jet vectored at an angle larger than the 
ramp deflection angle, the jet shear layer creates a second expansion region. 
Therefore, a second supersonic region may be present as a result of thrust-induced 
expansion. Figure 17 shows supersonic regions on this configuration. These 
regions are indicated by shading. As shown, a supersonic region occurs at the 
ramp even if the jet is not vectored. But the rest of the region about the nozzle 
remains subsonic. When the jet is ejected at an angle of 30 degree, not only the 
ramp region but also the region near the nozzle exit becomes supersonic. This 
calculation indicates that the flow outside the jet is deflected with the jet. 
Therefore, another expansion region is created by the thrust-vectoring jet. Figure 
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18 shows the Mach number distribution along the center line of the upper surface. 
For the case with the 30-degree jet, the Mach number distribution indicates that 
supersonic regions exist at the ramp and the jet exit plane. On the other hand, 
the Mach number distribution for the case with a straight jet shows only a 
supersonic region at the ramp. The rest of the region remains subsonic. As can 
be seen, at the region near the jet exit, the flow Mach number is even lower than 
the freestream mach number (M < 0.9 near the exit). Possible reasons for this 
phenomenon will be discussed in the next section. 

3.4 Jet Static-Pressure Effects 

Experimental work by Compton (Ref. 44) showed that the exhaust gas 
physical properties can be related to the jet plume shape and the entrainment and 
so influence the jet interference on afterbody pressure distribution. When the jet 
static pressure is greater than the outside pressure, the jet plume boundary will 
expand (Fig. 19). When the jet boundary expands, plume blockage occurs. Plume 
blockage can be related to the initial expansion angle of the plume boundary. The 
angle is the Prandtl-Meyer expansion angle of the jet flow when expanding from 
its internal exhaust static pressure to the external static pressure at the exit 
plane. As a result, a Prandtl-Meyer wedged-shaped region will be created at the 
nozzle edge, and creation of expansion waves is to be expected (Ref. 44). 

Expansion waves travel across the jet stream and reflect from the free jet 
boundaries as weak shocks. This wave-reflecting process goes on until viscous 
effects damp out the flow. With the waves reflecting inside the jet, alternating 
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high and low pressure regions will be created inside the jet. Two nozzles of NPR 
= 5 with different jet pressures are calculated to demonstrate the relation between 
the jet pressure and shock waves contained in the jet stream. The nozzle exit 
conditions are given in the following table. 
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Calculated velocity distribution for the first case is illustrated in Figure 20. The 
pressure and Mach number contours are presented in Figures 21 and 22, 
respectively. The first case, Mj et =1.3 and Pj e t/Pinf = L8 > P roduces strong shock 
waves inside the jet (Fig. 21a). The pressure contour on the symmetric plane 
clearly shows the alternation between low- and high-pressure regions. However, 
for the second jet condition (P jet = P inf , Fig. 21b), the pressure contours on the 
symmetric plane hardly show the shock waves pattern. Figure 22 shows the jet 
shear layers associated with these two cases. When the jet pressure is high, the 
jet plume tends to expand at the nozzle exit (Fig. 22a). When the plume blockage 
occurs, as a result of flow turning and deceleration, pressure distribution on the 
boattail will rise. As shown in Figure 23, the jet with a higher static pressure also 
induces higher pressure on the afterbody. These calculated results for straight 
jets are summarized in Figure 24. These calculations demonstrated the capability 
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of computing the effect of jet static pressure on flow properties of the shock- 
containing jet. 

Straight jet calculations show that a jet with high static pressure produces 
larger Prandtl-Meyer expansion and stronger jet-containing shock. As a result of 
higher jet pressure, the upper surface of the afterbody sustains higher pressure 
distribution. For the thrust vectoring application, when the jet is deflected 
downwards, low pressure on the upper surface is desired. However, for a deflected 
jet, if the upper portion of the jet plumes out, the upper portion of the jet stream 
does not follow the jet vectored axis. Therefore, the external upper surface flow 
decelerates and induces a more positive pressure region. To demonstrate the 
effect of the jet static pressure on the pressure distribution of the afterbody, two 
cases with the same NPR yet different static pressures are tested. Nozzle exit 
conditions of these cases are listed below. 
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The computed results are presented in Figures 25-27 for the Mach contours, 
pressure contours, and pressure distributions on the plane of symmetry, 
respectively. By comparing Figures 25a and 25b, it is seen that with a higher jet 
static pressure, flow expansion on the upper surface of the afterbody is much 
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reduced. This phenomenon can also be seen from the pressure contours shown in 
Figures 26a and 26b. Figure 27 again indicates that a more negative C p on the 
upper surface of the afterbody is obtained with a lower jet static pressure. 

To evaluate the aerodynamic performance of these configurations, computed 
total forces based on the wing reference area are compared as follows: 
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It is seen that the nozzle with a lower static jet pressure produces a slightly 
higher lift and more negative pitching moment. However, the drag coefficient is 
much higher. These calculations demonstrated static pressure effects on the 
vectored-jet performance. 

3.5 Delta Wing with Thrust Vectoring Effect 

This section discusses thrust vectoring effects on the leading edge vortices. 
Various theories of vortex formation have been proposed. Based on presently 
available results, a simple model (Ref. 46) for vortex formation will be used to 
explain the computational results. The leading edge vortex can be divided into 
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two parts: the inner core and the outer core. The inner core is characterized by a 
high axial velocity, which can reach three times the freestream speed. The outer 
core is characterized by axial and swirling velocity components that are 
approximately equal in magnitude. A model for the sequence of events leading to 
the formation of a leading edge vortex is described as follows. 

(1) The leading-edge sweep and pressure difference around the leading 
edge cause the flow to rotate about a central axis to form a rotating 
core. This core is called the inner core. 

(2) After the inner core is formed, the swirling component of the inner 
core causes the separated shear layer to wrap around the inner core 
and form the outer core. For conditions where the vortex breaks 
down, the inner core diffuses, resulting approximately equal 
magnitude in the swirling and axial velocity components. 

A 63-degree delta wing with a thick base (Fig. 28) is set up to investigate 
the vectored-jet effect on the leading-edge vortex. In this figure, shaded trapezoids 
indicate the nozzle exit planes. According to the sequence of leading-edge vortex 
formation, pressure difference and shear layer viscosity are two important factors. 
Previous delta wing calculations obtained by the CFL3D show good agreement 
with experimental data (Ref. 35). For the present configuration, the same grid 
topology as described in this reference is used. As shown in Figure 29, a block of 
60 x 21 x 25 points is constructed for the wing and wake region. A second block of 
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25 x 17 x 9 points is constructed to fill the space created by the thick base. A 
patched-grid method is used between these two main blocks. 

This configuration was computed under nozzle exit conditions listed below. 
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As shown in Figures 30-32, because of the thrust vectoring effect, the 
pressure distributions on the upper surface in all these cases is different in the 
trailing-edge regions. However, along the leading edge near the apex the pressure 
distributions are quite similar. Therefore, the vorticity near the apex for all cases 
is not much affected. As shown in Figure 33, inner core patterns of all cases are 
similar in the front portion of the vortex cores. However, near the trailing edge, 
the core with thrust vectoring effect tends to contract. One possible reason for 
this phenomenon is that the thrust momentum tends to bring the vortex core 
toward the jet. With the extra flow acceleration, the vortex core has more axial 
momentum along the streamwise direction. Since the outer core vortex is induced 
by the inner core swirling velocity, the outer vortex sheets should also be different. 
As shown in Figure 34, with higher momentum in the axial direction, the vortex 
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sheets in the case with vectored jet are more organized than the case without 
thrust vectoring effect. 

One way to present the vortex structure is by showing the sectional total 
pressure contours. For an irrotational flow, the total pressure is constant 
throughout the flow field. When a vortex flow occurs, the total pressure reduces 
as the vorticity increases. Because of the vorticity, the inner core bears a lower 
total pressure distribution compared with other regions. As shown in Figure 35, 
the vortex under the effect of a vectored jet has lower total pressure in the inner 
core than the configuration with a straight jet. In other words, a positively 
vectored jet (5 d > 0°) can increase the vorticity. As mentioned earlier, the vortex 
breaks down when the axial and swirling velocity components attain 
approximately equal magnitude. If the axial velocity decading rate can be 
reduced, the vortex breakdown can be delayed. The present computation shows 
that with the induced effect from a vectored jet, the axial velocity of a vortex core 
in the region near the nozzle can be increased. Thus, it is possible that a vectored 
jet may delay the occurrence of vortex breakdown. 
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4. CONCLUSIONS 


Aerodynamics induced by thrust vectoring schemes have been studied with 
a computational method. Specifically, configurations with 2D nozzles were used to 
investigate the aerodynamic effects related to thrust vectoring operations by thin- 
layer Reynolds-averaged Navier-Stokes solutions. To simulate the jet stream, jet 
velocity components and static pressures were specified at the nozzle exit. Two 
configurations were used in the calculation to verify the code, a rectangular jet in 
the crossflow, and a wing-body combination with a rectangular nozzle. Calculated 
results for these two configurations were compared with available data. 

Additional calculations were performed on a fuselage alone with a rectangular 
nozzle and a thick delta wing of 63-deg sweep to elucidate the thrust vectoring 
effects. The results indicated that: 

(1) Vectored thrust created a blockage effect in the flow field. A 
higher pressure region was created in the compression side, 
and a lower pressure was created in the expansion side. 

(2) A vectored jet induced flow expansions and shock waves in the 
jet at a transonic speed such that thrust vectoring effects were 
reduced. 

(3) A vectored jet induced thrust vectoring effects as shown in (1) 
on the nozzle as well as on near-by lifting surfaces. 

(4) A positively vectored jet made the swirling flow of a leading- 
edge vortex tighter, i.e., less diffusive. 
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Figure 1 Surface Velocity Vector Plot; 90-Degree Crossjet; R - 2.2 
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IL 

(a) Experiment (Taken from Ref. 22). 



(b) Computation 

Figure 3 Surface Pressure Contour Comparison; 90-Degree Crossjet; 
R = 2.2. 
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Figure 4 Geometry of the Generic Airplane Configuration. 
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(b) Locations of Pressure Orifices 



(c) Turning-Vane Geometry 


Figure 5 Nozzle and Turning-Vane Equipped in the Generic Airplane 
Configuration. 
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Figure 6 ' Grid Topologies of the Generic Airplane Calculations. 
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Figure 6 Continued. 
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Figure 7 


Effective Jet-Turning Angles (Taken from ref. 19). 
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Figure 8 I Jet Total Pressure Conversion from Thrust Coefficient 
(Taken from Ref. 19). 
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Figure 9 Pressure Contours on the Upper Surface with and without 
Jet Deflection. C T = 0.06, M„ = 0.8, a = 3.2 deg. 
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(b) Lower Surface 

Figure 11 Sectional Pressure Distributions; 2y/b = 0; M, 
3.2°; Mj el = 1.0; P jcl = P~. 









Figure 15 Geometry of a Fuselage with a Rectangular Nozzle. 


66 



Figure 16 Grid Topologies for the Fuselage Calculation. 
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(b) Second Block 


Figure 16 Continued. 
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(c) Third Block 


Figure 16 Concluded. 
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(b) 5 d = 30° 

Figure 17 Indication of Supersonic Regions on the Fuselage. 
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Figure 18 Sectional Mach Number Distributions on the Fuselage; 
a = 0°; M = 0.9; 2yft = 0.0; Mj et = 1.7; P jet = P inf . 
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Figure 19 Factors Influencing Jet Interference (Taken from Ref. 43). 
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Figure 20 Sectional Velocity Vectors of the Fuselage with a 

Rectangular Nozzle; = 1.3; P jct = 1.8 P inf ; M = 0.8; a = 0°; 
2y/b = 0. 
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Mach Number 


- M jrt = 1.7, = 1.0 P inf 



(a) Mach Number 



Figure 24 Flow Properties along the Centerline of the Jet, a — 0 , 
M = 0.8. 
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Figure 26 Pressure Contours of the Fuselage Configuration; 
a = 0°; 2y/b = 0. 
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Figure 26 Continued. 





Figure 28 Geometry of the 63-Degree Delta Body with Two 
Trapezoidal Nozzles. 
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(a) First Block 


(b) Second Block 



Figure 29 Grid Topologies for the Delta Body Calculation. 
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(a) Upper Surface 


Figure 30 


Pressure Distribution on the 63-Degree Delta Body without 
Jet; a = 25°; = 0.6. 
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(b) Lower surface 


Figure 30 Continued. 
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(b) Lower Surface 


Figure 31 Continued. 
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(b) Lower Surface 


Figure 32 Continued. 
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Figure 33 Inner Vortex Filament of the 63-Degree Delta Body with 
Various Jet Conditions. 



(a) Jet Off 



(b) 8(j « 0 deg. 



(c) 5d- 30 deg. 

Figure 34 Outer Vortex Filament of the 63-Degree Delta Body with 
Various Jet Conditions. 
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(a) 5d = 0 deg. 



(b) 5d = 30 deg. 


Figure 35 Total Cross Section Pressure Contours of the 63-Degree 
Delta Wing; a = 25°; M_ = 0.6; Mj Ct = 2; P jot = P inf . 
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APPENDIX A: 


GENERALIZED TRANSFORMATION 

Consider a completely general transformation of the form 


% - ^(x, y, z) 

T1 = ti(x, y, z) 

C = C(x, y, z) 


(A.1) 


The partial derivatives in the x, y, z coordinates can be written as 


3 _ 3 t 3 3 ^ 

3x 3^ 3r( 3£ 

3 3 t 3 3 i- 

T - - y + y + 3F^ y 
3y 3£ y 3r| * 3£ * 

3 _ 3 t 3 3 * 

3z 3£ 3n BC, 


(A.2) 


Also, the differential form of the generalized coordinate (£, T], Q can be written as 


d£ = £ x dx + ^ydy + E, Z dz 
drj = ri x dx + ti ydy + r\ z dz 

dC = Cx dx + C ydy + Cz dz 


(A.3) 


A.1 



In a matrix form, Equation (A.3) can be written as 




Sx £ y 

dx 

drj 

= 

Tlx Tly Tlz 

dy 

LdeJ 


C* C y 

dz 


(A.4) 


In a similar manner, the differential form of (x, y, z) can be written as 


dx 


X E, x n x^ 

d*f 

dy 

= 

y^ yq y^ 

drj 

_dz_ 



dC. 



Z C. 



(A.5) 


Therefore, 


4x 5y 5* 


* 

X 

►f 

WLP 

, X 

"H X Tf y TJ z 

= 

y^ yq y; 

?x Cy tz_ 


Zq Z C_ 


= J 


(y n z ; 

-( y 

. ( ^ z n 


- y^) 

- y^) 

- yqZ{0 


-(Xq z c - x c z n ) 

(x^Z£ - x^z^) 

-(x^Zq - Xq 2 ^ 


(Xqy; - x^y n ) 
-(x^yxub^ - x^y^) 
(x^yq - x n y^) 


(A.6) 


A.2 



j = [(yn z c - (y^ z c - yc z ^> x ri + x c ( y^ z ri - yn z ^l 


Equation (2) can be written as 

-> 03. 
at 

♦ 4 (F - F u«x ♦ 4 (F - F u>nx + 4 (F - F U K X 

a^ an a; 

*JL(G - G v ) * 4r (G - G u>ly + 4 (G - G uKy 
dr i or | 

♦ 4< H - H u«x + 4-® - H v>ny * 4® - H vKz = 0 

dc, an ac, 


where 


aQ_aQa^ + aQan + aQa^ + aQ 
at a^ at an at ac at at 


Equations (8) and (9) can be obtained by rearranging Equation (A.7). 


A.3 



APPENDIX B: 


DERIVATION OF VAN LEER’S FLUX VECTOR SPLITTING 

FORMULATION 

The goal of van Leer’s flux vector splitting (Ref. 37) is to split the flux ffBJ) 
into a forward flux f’XtD) and a backward flux ffGJ); that is, 
flco) = f + (oo) + f"(co); and 

(2) df + /dco have all non-negative eigenvalues, and 


df“/dco have all non-positive eigenvalues. 


Considering the one-dimensional Euler equations, the full flux can be 
written as 


f(p, u, p, e) 


( \ 


pu 

pu 2 + p 
j)ue + pu ; 


= flp, C, M) = 


fpCM 


pC 2 (M 2 + 

pC 3 M(lM 

2 



A 



J 


The flux is split under the following restriction: 


B.l 



( 1 ) 


f*(co) = flco) when M > 1 
f-(co) = flco) when M < -1 


(2) f*(M) = ± fl-M), if KM) = fl-M) 


(3) 


df* 

dco 


must be continuous 


(4) 


df ± 

dw 


must have one eigenvalue vanish for |M j < 1 


(5) f^M) must be a lowest-possible-degree polynomial in M. 


The polynomials of the flux are determined as 



1 


± pC(M ± l) 2 /4 

f* = 


= 

£f[(y - l)u ± 2C]/y 


y 


ff[(Y - l)u ± 2C] 2 /2(y 2 - 1) 


B.2 



For three-dimensional Euler equations, the x-flux can be written as 


f = 


■s 


pu 


pu^+p 


puv 

— 

puo) 


^(exp)u j 

V 


pu 

p(u 2 + C 2 /y) 

puv 

puo) 

pu[JL(u 2 + v 2 + w 2 ) + C 2 /y-l] 
2 


where 


e = 


y - 1 2 


+ -lp(u 2 + v 2 + w 2 ), C = 


N 


yp 

P 


f = 



“ 

ft 


f? 


ft 

= 

ft 


ft 



± pC(M x ± l) 2 /4 
ff(y -l)u ± 2C)/y 


f l ‘ v 

ff • 0) 


' 

r 

2 


(y - l)u ± 2C 

1 , 2 2D 

+ — (v + w ) 

2 


1 

to 

to 

1 

M 


B.3 



In Equation (25), 


F = Mi [$*F * i y G + ? Z H] 

J 


F± = Mi [^F* + 5 G± * $ 2 H±] 

J J 


±pC(M^ ± l) 2 /4 

f^(u| x + v^ y + w£ z ) + U(f?[-u * 2C]/y) 
f^(uf| x + vf| y + wf) z ) + ^ y (f^[-u ± 2C]/y) 
fj(u£x + v£y + w£ z ) + $,(ff[-u ± 2C]/y ) 
fj{[-(y - l)u 2 ± 2(y - l)uC + 2a 2 ]/(y 2 + 1) 
+ (u 2 + v 2 + w 2 )/2} 



B.4 



APPENDIX C: 


BOUNDARY CONDITIONS USED IN CFL3D 
Standard Boundary Conditions: 


1000 

free stream 

1001 

reflection of x-z plane 

1002 

extrapolation 

1003 

inflow/outflow 

1004 

viscous surface 

1005 

inviscid surface 

1006 

x-y plane 

1007 

axi symmetric 


Special Purpose Conditions: 
1146 


jet exit plane 


C.l 




